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ABSTRACT 

Future weak lensing surveys will directly probe the density fluctuation in the universe. 
Recent studies have shown how the statistics of the weak lensing convergence field is 
related to the statistics of collapsed objects. Extending earlier analytical results on 
the probability distribution function of the convergence field we show that the bias 
associated with the convergence field can directly be related to the bias associated with 
the statistics of underlying over-dense objects. This will provide us a direct method 
to study the gravity induced bias in galaxy clustering. Based on our analytical results 
which use the hierarchical ansatz for non-linear clustering, we study how such a bias 
depends on the smoothing angle and the source red-shift. We compare our analytical 
results against ray tracing experiments through N-body simulations of four different 
realistic cosmological scenarios and found a very good match. Our study shows that the 
bias in the convergence map strongly depends on the background geometry and hence 
can help us in distinguishing different cosmological models in addition to improving 
our understanding of the gravity induced bias in galaxy clustering. 

Key words: Cosmology: theory -- large-scale structure of the Universe - Methods: 
analytical 



1 INTRODUCTION 

I Weak gravitational lensing (Bartelmann & Schneider 1999) is the distortion of images of the background galaxies through the 
tidal gravitational field of large scale background mass distribution. Study of such distortions provides us an unique way to 
probe statistical properties of intervening mass distribution. 

Although there is no generally applicable definition of weak lensing the common aspect of all studies in weak gravitational 
lensing is that measurements of its effects are statistical in nature. Pioneering works in this direction were done by Blandford 
et aL (1991), Miralda-Escude (1991) and Kaiser (1992) based on the earlier work by Gunn (1967). Current progress in weak 
lensing can broadly be divided into two distinct categories. Where as Villumsen (1996), Stebbins (1996), Bernardeau et al. 
(1997) and Kaiser (1987) have focussed on linear and quasi-linear regime by assuming a large smoothing angle, several authors 
have developed a numerical technique to simulate weak lensing catalogs. Numerical simulations of weak lensing typically 
employ N-body simulations, through which ray tracing experiments are conducted (Schneider & Weiss 1988; Jarosszn'ski et 
al. 1990; Lee & Paczyn'ski 1990; Jarosszn'ski 1991; Babul & Lee 1991; Bartelmann & Schneider 1991, Blandford et al. 1991). 
Building on the earlier work of Wambsganns et al. (1995, 1997, 1998) most detailed numerical study of lensing was done by 
Wambsganns, Cen & Ostriker (1998). Other recent studies using ray tracing experiments have been conducted by Premadi, 
Martel & Matzner (1998), van Waerbeke, Bernardeau & Mellier (1998), Bartelmann et al (1998) and Couchman, Barber & 
Thomas (1998). While peturbative analysis can provide valuable information at large smoothing angle such analysis can not 
be used to study lensing at small angular scales as the whole perturbative series starts to diverge. 

A complete analysis of the weak lensing statistics at small angular scales will require a similar analysis for the underlying 
dark matter distribution which we do not have at present. However there are several non-linear ansatze which predict a tree 
hierarchy for matter correlation functions and are thought to be successful to some degree in modeling results from numerical 
simulations. Most of these ansatze assumes a tree hierarchy for higher order correlation functions and they disagree with each 
other by the way they assign weights to trees of same order but of different topologies (Balian & Schaeffer 1989, Bernardeau 
& Schaeffer 1992; Szapudi & Szalay 1993). Evolution of two-point correlation functions in all such approximations are left 
arbitrary. However recent studies by several authors (Hamilton et al 1991; Nityananda & Padmanabhan 1994; Padmanabhan 
et al. 1996; Jain, Mo & White 1995; Peacock & Dodds 1996) have provided us a very accurate fitting formula for the evolution 
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of the two-point correlation function which can be used in combination with these hierarchical ansatze to predict clustering 
properties of dark matter distribution in the universe. 

Most recent studies in weak lensing have mainly focussed on lower order cumulants (van Waerbeke, Bernardeau & 
Mellier 1998, Hui 1999, Munshi & Jain 1999a), cumulant correlators (Munshi & Jain 1999b) and errors associated with their 
measurement from observational data using different filter functions (Roblinsky et al. 1999). However it is well known that 
higher order moments and more sensitive to the tail of the distribution function which they represent and are sensitive to 
measurement errors due to finite size of catalogs (Szapudi & Colombi 1996). On the other hand numerical simulations involving 
ray tracing techniques have already shown that while probability distribution function (PDF) associated with density field is 
not sensitive to cosmological parameters its weak lensing counterpart can lift such a degeneracy and can help us in estimation of 
cosniological parameters from observational data (Jain, Scljak & White 1999). Munshi & Jain (1999a, b) extended such studies 
to show that the hierarchical ansatz can actually be used to make concrete analytical predictions for lower order statistical 
properties of the convergence field. Valageas (1999a) has used hierarchical ansatz for computing error involved in estimation 
of fJo and A from SNcIa observations. A similar fitting function was recently proposed by Wang (1999). Our formalism is 
similar to that of Valageas (1999a) although the results obtained by us include the efi^ect of smoothing and analysis presented 
here correspond to the joint probability distribution function of two smoothed patches and hence the bias associated with 
them. Our formalism presented here is an extension of the study of the PDF using similar techniques (Munshi & Jain 1999b, 
Valageas 1999b), it can very easily be extended to multi-point PDF and hence to compute the Sn parameters associated 
with "hot spots" of the convergence maps. Analytical results and detailed comparison against numerical simulations will be 
presented elsewhere (Munshi & Coles 1999a,b). 

The paper is organized as follows. In section 2 we describe briefly ray tracing simulations. Section 3 presents most of the 
analytical results necessary for computation of the bias of the smoothed convergence field Ks. Results of data analysis are 
presented in section 4. Section 5 is left for discussion of our result in general cosmological framework. 



2 GENERATION OF CONVERGENCE MAPS FROM N-BODY SIMULATIONS 

Convergence maps are generated by solving the geodesic equations for the propagation of light rays through N-body simulations 
of dark matter clustering. The simulations used for our study are adaptive P^M simulations with 256^ particles and were 
carried out using codes kindly made available by the VIRGO consortium. These simulations can resolve structures larger than 

30h~^kpc at z = accurately. These siunilations were carried out using 128 or 256 processors on CRAY T3D machines at 
Edinburgh Parallel Computer Center and at the Garching Computer Center of the Max-Planck Society. These simulations 
were intended primarily for studies of the formation and clustering of galaxies (KaufFmann et al 1999a, 1999b; Diaferio et 
al 1999) but were made available by these authors and by the Virgo Consortium for this and earlier studies of gravitational 
lensing (Jain, Scljak & White 1999, Reblirisky et al. 1999, Munshi & Jain 1999a,b). 

Ray tracing simulations were carried out by Jain et al. (1999) using a multiple lens-plane calculation which implements 
the discrete version of recursion relations for mapping the photon position and the Jacobian matrix (Schneider & Weiss 1988; 
Schneider, Ehler & Falco 1992). In a typical experiment 4 x 10® rays are used to trace the underlying mass distribution. The 
dark matter distribution between the source and the observer is projected onto 20 - 30 planes. The particle positions on each 
plane are interpolated onto a 2048'^ grid. On each plane the shear matrix is computed on this grid from the projected density 
by using Fourier space relations between the two. The photons are propagated starting from a rectangular grid on the first lens 
plane. The regular grid of photon position gets distorted along the line of sight. To ensure that all photons reach the observer, 
the ray tracing experiments are generally done backward in time from the observer to the source plane at red-shift z = Zs- 
The resolution of the convergence maps is limited by both the resolution scale associated with numerical simulations and also 
due to the finite resolution of the grid to propagate photons. The outcome of these simulations are shear and convergence 
maps on a two dimensional grid. Depending on the background cosmology the two dimensional box represents a few degree 
scale patch on the sky. For more details on the generation of /t-maps, see Jain et al (1999). 



3 ANALYTICAL PREDICTIONS 

In this section we will provide necessary theoretical background for the bias function in the context of hierarchical clustering. 
We will be using the following line element for the background geometry of the universe: 

dr^ = -cdt'^ + a^{t){dx' + r^(x)d^n) (1) 

Where we have denoted angular diameter distance by r{x) and scale factor of the universe by a{t). r{x) = K^^^^ sin(if ~^''^x) 
for positive curvature, r(x) = {--K)~^^^ s\a\\{{—K)~^^^x) for negative curvature and x for the flat universe. For a present 
value of of Ho and f2o we have K = (f2o — ^)Ho. The various parameters charecterising different cosmological models are 
listed in Table - 1. 



Table 1. Cosmological Pcirameters charecterising different models 





SCDM 


TCDM 


LCDM 


OCDM 


r 


0.5 


0.21 


0.21 


0.21 


no 


1.0 


1.0 


0.3 


0.3 


Ao 


0.0 


0.0 


0.7 


0.0 




0.6 


0.6 


0.9 


0.85 


^^0 


50 


50 


70 


70 



3.1 The Formalism 

The statistics of weak lensing convergence k is very much similar to that of the projected density field. In what follows we will 

be considering a small patch of the sky where we can use the plane parallel approximation or small angle approximation to 
replace spherical harmonics by Fourier modes. The 3D density contrast 5 along the line of sight when projected into 2D sky 
with a weight function ui(x) will provide us the projected density contrast or the weak- lensing convergence at a direction 7. 



«(7i) = / dxMxi)S{r{xhi) (2) 
Jo 

In all our discussion we will be placing the sources at a fixed red-shift (an approximation not too difficult to modify for more 
realistic description), the weight function can be expressed as a;(x) = S/4ac~'^ Hoflmr{x)r{xs — x)/^(Xs)- Where Xs is the 
comoving radial distance to the source at a redshift Zg- Fourier decomposition of 6 can be written as: 

(3) 



Where we have used fey and ki_ to denote components of wave vector k parallel and perpendicular to line of sight direction 7. 

In small angle approximation however one assumes that k± much larger compared to fcy . We will denote the angle between 
the line of sight direction 7 and the wave vector k by 9. Using the definitions we have introduced above we can compute the 
smoothed projected two-point correlation function (Peebles 1980, Kaiser 1992, Kaiser 1998): 

(^(7i)«(72))c = ^"dxi^ j eMOl) P{^^)wi{Wo). (4) 

Where we have introduced a new notation I = r(x)kx which denotes a scaled wave vector projected on the surface of the sky. 
The average of two-point correlation function (ks) smoothed over an angle 9o with a top-hat smoothing window W2{lOo) is 
useful to quantify the fluctuations in Ke which is often used to reconstruct the matter power spectrum i-'(k) (Jain, Selzak & 
White 1998). 

Similar analysis for the higher order cumulant correlators (Szapudi & Szalay 1997, Munshi & Coles 1999) of the smoothed 
convergence field relating (Ky^(7i)«;"(72))c with multi-spectra of underlying dark matter distribution Bp (Munshi & Coles 
1999a): 

(K?(7i)-.(72))c = p ^dx J ^ j ^,W,{heo)W,{l,eo)W,{heo)eMihe^,)B, , (6) 



(k^(7i)ks(72))c 
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(«?(7i)«?(72))c = p ^'^X j ^ J ^ J ^W{heo)W2{l2eo)W2{l30o)W2{l49o)eMiih+h)Oi2) 

bJ-^ (8) 
* \r{x)' r{x)' r{x)' r{x) J ' 

In general we can express the cumulant correlators of arbitrary order (Kr(7i)t?(72))c in terms of multi-spectra Bm+n as: 



We will use and extend these results in this paper to show that it is possible to compute the whole bias function 6(> Ks), 
i.e. the bias associated with those spots in convergence map which Ks is above certain threshold (which acts as a generating 
function for these cumulant correlators) from the statistics of underlying over-dense dark objects. Details of analytical results 
presented here can be found in Munshi & Coles (1999b). 

3.2 Hierarchical Ansatze 

In deriving the above expressions we have not used any specific form for the matter correlation hierarchy, however the length 
scales involved in small angles are in the highly non-linear regime. Assuming a tree model for the matter correlation hierarchy 
in the highly non-linear regime one can write the most general case as (Groth & Peebles 1977; Pry & Peebles 1978; Davis & 
Peebles 1983; Bernardeau & Schaeffer 1992; Szapudi & Szalay 1993): 

(N-l) 

U(ri,...rN)= Qn,c J2 n ^(''''''J) 

a, N— trees labellings edges(i,j) 

It is interesting to note that an exactly similar hierarchy develops in the quasi-lineax regime in the limit of vanishing 
variance (Bernardeau 1992), however the hierarchal amplitudes Qjv.a become shape dependent in such a case. In the highly 
nonlinear regime there arc some indications that these functions become independent of shape parameters as has been proved 
by studies of lowest order parameter Qs = Q using high resolution numerical simulations (Sccociamarro et al. 1998). In the 
Fourier space such an ansatz will mean that the whole hierarchy of multi-spectra Bp can be written in terms of sum of products 
of power-spectra, e.g. in low orders we can write: 

B2(kl, k2, ka)^ = Q(P(kl)P(k2) + P(k2)P(k3) + P(k3)P(ki)), (11) 

P3(ki, k2, ka, k4)^ = P„P(ki)P(ki + k2)P(ki -|- ka + k3) + cyc.perm. -|- i?(,P(ki)P(k2)p(k3) + cyc.perm. (12) 

Different hierarchal models differ the way they predict the amplitudes of different tree topologies. Bernardeau & Schaeffer 
(1992) considered the case where amplitudes in general are factorizable, at each order one has a new "star" amplitude and 
higher order "snake" and "hybrid" amplitudes are constructed from lower order "star" amplitudes (see Munshi, Melott & 
Coles 1999a,b,c for a detailed description). In models proposed by Szapudi & Szalay (1993) it is assumed that all hierarchal 
amplitudes of a given order are actually degenerate. 

We do not use any of these specific models for clustering and only assume a hierarchal nature for the higher order 
correlation functions. Galaxy surveys have been used to study these ansatze. Our main motivation here is to show that 
weak-lensing surveys can also provide valuable information in this direction, in addition to constraining matter power-spectra 
and the background geometry of the universe. The most general form for the lower order cumulant correlators in the large 
separation limit can be expressed as: 



{Ks{'yi)Ksiy2)}c = 2(3aCa[«;9o«;9i2] = C'2iC3[«;9oKei2] = CJi(/ti,)c{«;s(7i)«s(72))c, (13) 

(«s(7i)ks(72))c = (3i?a -|-6i?6)C4[Keo/tei2] = C^iC4[KeoKei2] = C3i(/ts)?(«s(7i)«s(72))c, (14) 

(Ks(7i)k?(72))c = 4RbC4[KlgKei2] = ^22^41^001^812] = C22(l4)c(l^si'yi)Ks{-y2))c, (15) 

(K^(7i)«:,(72)),: = {24Sa + 36St + 4S.)C,[4o^<'e,2] = C2^C,[kI,ko,2] = Cri{K')?{AC,(7l)K,(72))c, (16) 

(«?(7i)'^'(72))c = (125„+6S6)C5[4„«ei2] = C73V5[«eo«ei2] = C:52(^')c('^.(7i)'^.(72))c. (17) 

Where denotes the cumulant correlators of the convergence field and denotes the cumulant correlators for the 
underlying mass distribution. Extending above results to arbitrary order we can write: 

(«?(7i)«I(72))c = C;,Cp+,[K<f+''-^\e,2] = C;,(«?)i^+^-''(K.(7i)«.(72))c. (18) 
where C^q denotes the cumulant correlators for the underlying mass distribution, 

Ct[K^Kei2]=^ ^.2^-if|^) ^fli^^gi2c'x. (19) 



similarly the following notations were used to simplify the above expressions: 



f^p(- 



««o = / 7^P{ -^)whieo), (20) 

"'^^ ^ j ^^^^ 

The hierarchical expression for the lowest order cumulant i.e. ^3 was derived by Hui 1998. He also showed that his result 
agrees well with numerical ray tracing experiments by Jain, Seljak and White (1998). More recent studies have shown that 
higher order cumulants and even the two-point statistics such as cumulant correlators can also be reliably modeled in a similar 
way (Munshi & Coles 1999, Munshi & Jain 1999a). We extend such results in this paper to compute the complete bigis function 
&«(ks) (which is related to the low order two-point statistics such as the cumulant correlators defined above) associated with 
high /ts spots in the convergence map, 

Pk(«;i, K2)dKidK2 = Pk(ki)p,^(k;2)(1 + hn{iii)it2bK{n2))dK-i_dn2, (22) 

and its relation to the bias associated with collapsed objects in underlying density field 1 + 5. 

3.3 The bias associated with collapsed objects 

The success of analytical results in predicting the lower order cumulants, cumulant correlators and the one-point smoothed 
PDF(Munshi & Coles 199a,b; Munshi & Jain 1999a,b; Valageas 1999a,b), motivates a general analysis of the bias associated 
with the high Ks spots in convergence maps and their relation with bias associated with the high peaks in the underlying 
mass distribution. For this purpose we found that formalism developed by Balian & Schaeffer (1989) and Bernardeau & 
Schaeffer (1992) (later extended Bernardeau (1992, 1994)) to be most suitable. These results are based on a very general tree 
hierarchy for higher order correlation function and the assumption that amplitudes associated with different tree-topologies 
are constant once in the highly non-linear regime. Later these results were generalized by Bernardeau (1992, 1994) for the case 
of quasi-linear regime too, where pcrturbativc dynamics can be used to make more concrete predictions. Errors associated 
with top-hat smoothing and its use with hierarchical ansatz has been studied in detail by Szapudi et al.(1992) and Boschan et 
al. (1994) and will be assumed small even in the case of two-point cumulant correlators throughout this paper. In this section 
we review the basic results from scaling models in the highly non-linear regime and quasi-linear regime before extending 
such models to the statistics of smoothed convergence field k{6o) where we will show that although weak lensing statistics 
probes the highly non-linear regime, projection effects make the variance smaller than unity justifying the use of quasi-linear 
results even though the generating function remains the one for highly non-linear regime. We will be using the small angle 
approximation in our derivation. Our results can in principle be generalized for the case of projected galaxy catalogs too, a 
detailed analysis will be presented elsewhere. 

3.3.1 The Generating Function 

In scaling analysis of the probability distribution function (PDF) the void probability distribution function (VPF) plays most 
fundamental role, which can related to the generating function of the cumulants or Sn parameters, (j){y) (White 1979, Balian 

& Schaeffer 1989) 

P.(0) = exp(-JVa(iVc)) = exp(-ffi^). (23) 

Where (0) is the probability of having no "particles" in a cell of of volume v, N is the average occupancy of these "particles" 
and Nc = N^2. The VPF is meaningful only for discrete distribution of particles and can not be defined for smooth density 
fields such as 5 or k{0o). However the scaling functions defined above a{y) = — are very much useful even for continuous 
distributions where they can be used as a generating function of one-point cumulants or Sp parameters, 

00 

-^(f ) = E ^y"- (24) 

The function (j>{y) satisfies the constraint Si = S2 = 1 necessary for proper normalization of PDF. The other generating 
function which plays a very important role in such analysis is the generating function for vertex amplitudes Un, associated 
with nodes appearing in tree representation of higher order correlation hierarchy (Qs = 1/2, Ra = and Rb = 1^3). 

G{r) = l-r+^r^-^r^ + ... (25) 

A more specific model for G{t) can be used, which is useful to make more specific predictions (Bernardeau & Schaeffer 1979): 



e(r)=(l + f) ^\ (26) 

Wo will relate Ka with other parameters of scaling models. While the definition of VPF do not use any specific form of 
hierarchical ansatz it is to realize that writing the tree amplitudes in terms of the weights associated with nodes is only 
possible when one assumes a factorizable model for tree hierarchy (Bernardeau &: Schaeffer 1992) and other possibilities which 
do not violate the tree models are indeed possible too (Bernardeau & Schaeffer 1999). The generating functions for tree nodes 
can be related to the VPF by solving a pair of implicit equations (Balian & Schaeffer 1989), 

m=vQ{r)-\yr-^G{r), (27) 

r = -y^Gir). (28) 

However a more detailed analysis is needed to include the effect of correlation between two or more correlated volume 

element which will provide information about bias, cumulants and cumulant correlators of these collapsed object (as opposed 
to the cumulants and cunmlant correlators of the whole convergence map, e.g. Munshi & Jain (1999a)). However we will only 
quote results useful for measurement of bias from ray-tracing simulations as detailed derivations of related results including 
related error analysis can be found elsewhere (Bernardeau & Schaeffer 1992, 1999; Munshi et al. 1999a,b,c; Coles et al. 1999). 

Notice that r(j/) (also denoted by f3{y) in the literature) plays the role of generating function for factorized cumulant 
correlators Cpi {Cpq = CpiCqi): 

-(^) = E (29) 

p=i 



3.3.2 The Highly Non-linear Regime 

The PDF p{5) and bias b{S) can be related to their generating functions VPF 4>{y) and t(j/) respectively by following equations 
(Balian & Schaeffer 1989, Bernardeau & Schaeffer 1992, Bernardeau & Schaeffer 1999), 



p{5) 
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(1 + 6)y - 4>{y) 



b{5)p{5) 



dy , . 
^r(y)exp 



(l + %-0(:y) 



6 



(30) 



(31) 



ft is clear that the function 0(j/) completely determines the behavior of the PDF p(5) for all values of 5. However different 
asymptotic expressions of c^(j/) govern the behavior of for different intervals of 5. For large y we can express (^{y) as: 



0(2/) = oy 



(32) 



Where we have introduced a new parameter w for the description of VPF. This parameter plays a very important role 
in scaling analysis. No theoretical analysis has been done so far to link lo with initial power spectral index n. Numerical 
simulations are generally used to fix w for a specific initial condition. Such studies have confirmed that the increase in power 
on smaller scales increases the value of ixi. Typically initial power spectrum with spectral index n = —2 (which should model 
CDM like spectra we considered in our simulations at small length scales) produces a value of .3 which we will be using in 
our analysis of PDF of the convergence field Ks (Colombi et. al. (1992, 1994, 1995). The VPF ^(y) and its two-point analog 
r(xj) both exhibit singularity for small but negative value of j/s. 



■r{y) 



- asr{u)s){y - ys) 

b,{y-ys)--^-\ 



(33) 



For the factorizable model of the hierarchical clustering the parameter lOs takes the value —3/2 and a, and 6s can be expressed 
in terms of the nature of the generating function G{t) and its derivatives near the singularity Ts (Bernardeau & Schaeffer 
1992): 



r(-i/2) 

2g'(r,)g"(r. 



g'{rsW"{n) 

1/2 



6"'{t,) 



2g'{Ts)g"{rs) 

e"'(rs) 



3/2 



(34) 
(35) 



As mentioned before the parameter ka which we have introduced in the definition of ^(t) can be related to the parameters 
a and uj appearing in the asymptotic expressions of ^(y) (Balian & Schaeffer 1989, Bernardeau & Schaeffer 1992), 



iU = ka/{ka+2), (36) 

a = t^k'^^/>^^+\ (37) 

Similarly the parameter ys which describe the behavior of the function (j}{y) near its singularity can be related to the 
behavior of Q{t) near Ts which is the solution of the equation (Balian & Schaeffer 1989, Bernardeau & Schaeffer 1992), 

''Q"{r^y ^^^^ 
finally we can relate ka to ys by following expression (see eq. (^)): 

= ^rTT-V (39) 
Q'{rs) 

or 

l-^-M^. + 2)-- ^^^^ 
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The newly introduced variable Xi, will be useful to define large 5 tail of the PDF p[5) and the bias b{5). Different asymptotes 
in 0(2/) are linked with behavior of p{5) for various regimes of 5. For very large values of variance i.e. ^2 it is possible to 
define a scaling function p{5) = 4ih{x) which will encode the scaling behavior of PDF, where plays the role of the scaling 

variable and is defined as ^j-^- We list different ranges of 5 and specify the behavior of p(5) and 6(5) in these regimes (Balian 



Schaeffer 1989). 



^2 

(l-")/2 



.,,.|(l|^)»p(-ii^). H.^-^lii^ 

The integral constraints satisfied by scaling function are Si — xh{x)dx = 1 and S2 = x^h{x)dx — 1. These take 
care of normalization of the function p{S). Similarly the normalization constraint over b{S) can be expressed as Cii = 
JJ" xh{x)h{x)dx — 1, which translates into d5b{5)p{5) = and dS5b{5)p{5) — 1. Several numerical studies have been 
conducted to study the behavior of h{x) and b{x) for different initial conditions (e.g. Colombi et al. 1992,1994,1995; Munshi 
et al. 1999, Valageas et al. 1999). For very small values of S the behavior of p{5) is determined by the asymptotic behavior of 
(jiiy) for large values of y, and it is possible to define another scaling function g{z) which is completely determined by ui, the 
scaling parameter can be expressed as z = (1 + 6)0'^^^^""^^ ^2^^^ However numerically it is much easier to determine uj 
from the study of fj{y) compared to the study of g{z) (e.g. Bouchet & Hernquist 1992). 



l + 5«6; piS)^a^--^^f-JJl-;^e^p 



To summarize, we can say that the entire behavior of the PDF P{S) is encoded in two different scaling functions, h(x) and 
g(z) and one can also study the scaling properties of b{S) in terms of the scaling variables x and 2; in a very similar way. These 
scaling functions are relevant for small and large 5 behavior of the function p{S) and b{5). Typically the PDF p{d) shows a 
cutoff at both large and small values of S and it exhibits a power-law in the middle. The power law behavior is observed when 
both g{z) and h{x) overlap and is typical of highly non-linear regime. With the decrease in ^2 the range of 5 for which p{5) 
shows such a power law behavior decreases finally to vanish for the case of very small variance i.e. in the quasi-linear regime. 
Similarly the bias is very small and a slowly varying function for moderately over dense objects but increases rapidly for over 
dense objects which are in qualitative agreement with PS formalism. 



3.3.3 The Quasi-linear Regime 

In the quasi-linear regime, a similar formalism can be used to study the PDF. However the generating function now can be 
explicitly evaluated by using tree-level perturbative dynamics (Bernardeau 1992; Bernardeau 1994). It is also possible to take 
smoothing corrections into account in which case one can have explicit expression of lj in terms of the initial power spectral 
index n. In general the parameters ka or lo charecterising VPF or CPDF are different from there highly non-linear values. 

For the purpose of weak lensing calculations it is important to notice that although the generating function for the matter 
correlation hierarchy for very small angular smoothing is the one from the highly non-linear regime, the analytical results 
that are useful are the ones from the quasi-linear regime as the variance of projected density field remains very small even for 
small smoothing scales. 

The PDF and bias now can be expressed in terms of Gs{t) (Bernardeau 1992, Bernardeau 1994): 





V 
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Figure 1. The cumulative bias 6r)(> rj) assoeiated with the reduced convergence field r; is plotted as a function of r) for three different 
values of variance (bottom to top) i.e. -^/i^ = .25, .5, and 1. Scaling parameter oj = 0.3 is used through out out our analysis. 
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Gsir) = G{t) - 1 ^ 5. (45) 

The above expression is valid for 5 < 5c where the 5c is the value of 5 which cancels the numerator of the pre-factor of the 
exponential function appearing in the above expression. For 5 > Sc the PDF develops an exponential tail which is related to 
the presence of singularity in <f){y) in a very similar way as in the case of its highly non-linear counterpart (Bernardeau 1992, 

Bernardeau 1994). 
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In recent studies it was shown that the underlying smoothed density PDF and its convergence counterpart Pk(ks) are 
exactly same (Munshi & Jain 1999b) under small angle approximation. In this paper we will extend such result to show that 
even the bias associated with the convergence map can be related to the bias associated with over-dense objects in a very 
similar manner. Such results can also be obtained from the extended PS formalism and we plan to present details of such 

analytical results elsewhere. 

It may be noted that similar analytical expressions for the PDF and bias can also be derived for the case of approximate 
dynamics sometime used to simulate gravitational clustering in the weakly non-linear regime or to reconstruct the projected 
density maps from convergence maps with large smoothing angles (e.g. Lagrangian perturbation theory which is an extension 
of Zeldovich approximation (Munshi et al. 1994)). 

3.4 The Bias of the Convergence Field 

To compute the bias associated with the peaks in the convergence field we have to fist develop an analytic expression for the 
generating field I3{yi,y2) for the convergence field Ks- For that we will use the usual defirntion for the two-point cumulant 
correlator Cpq for the convergence field (for a complete treatment of statistical properties of Ks see Munshi & Coles, 1999b). 

"^^^ - (fef>f+«-2(K.(7l)K.(72)> ~ ^ ' 

We will show that like its density field counterpart the two-point generating function for the convergence field Ks can also be 
expressed (under certain simplifying assumptions) as a product of two one-point generating functions P{y) which can then be 
directly related to the bias associated with "hot-spots" in the convergence field. 

/3.(j/i,J/2) = ^ ^ylyl = J2 E ^2^2 = Myi)My2) = r42/i)r4j/2) (48) 
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Figure 2. The construction of the bias function 6(/ts) for the smoothed convergence field from the bias 6(??s) associated with the 
smoothed reduced convergence field rjs. 



It is clear that the factorization of generating function actually depend on the factorization property of the cumulant correlators 
i.e. C^g — CpiC^j^. Note that such a factorization is possible when the correlation of two patches in the directions gammai 
and 72 {ns{"fi)K-s{"t2))c is smaller compared to the variance (/tg) for the smoothed patches 
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We will now use the integral expression for cumulant correlators (Munshi & Coles 1999a) to express the generating function 
which in turn uses the hierarchical ansatz and the far field approximation as explained above 
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It is possible to further simplify the above expression by separating the summation over dummy variables p and q, which will 
be useful to establish the factorization property of two-point generating function for bias Piyi,y2)- 
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We can now decompose the double sum over the two indices into two separate sums over individual indices. Finally using the 
definition of the one-point generating function for the cumulant correlators we can write: 
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The above expression is quite general and depends only on the small angle approximation and the large separation approxima- 
tion and is valid for any given specific model for the generating function G(r). However it is easy to notice that the projection 
effects as encoded in the line of sight integration do not allow us to write down the two-point generating function f3i^{yi,y2) 
simply as a product of two one-point generating functions I3n{y) as was the case for the density field 1 + 5. 

As in the case of the derivation of the probability distribution function for the smoothed convergence field Hs it will be 
much easier if we define a reduced smoothed convergence field r?s . The statistical properties of rjs are very similar to that of 
the underlying 3D density field (under certain simplifying approximation) and are roughly independent of the background 
geometry and dynamics of the universe, 



1 + 



(53) 



Where the minimum vale of i.e. k-min is defined as: 
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It is easy to notice that the minimum value of the convergence field will occur in those line of sight which are completely 
devoid of any matter i.e. 6 = —1 all along the line of sight. We will also find out later that the cosmological dependence of the 
statistics of Ks field is encoded in kmin and this choice of the new variable rjs makes its related statistics almost independent 
of the background cosmology. Repeating the above analysis again for the r), field, we can express the cumulant correlator 
generating function for the reduced convergence field rjs as: 
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While the above expression is indeed very accurate and relates the generating function of the density field with that of the 
convergence field, it is difficult to handol for any practical purpose. Also it is important to notice that the scaling functions 
such as h{x) for the density probability distribution function and h(x) for the bias associated with over-dense objects are 
typically estimated from numerical simulations specially in the highly non-linear regime. Such estimations are plagued with 
several uncertainties such as finite size of the simulation box. It was noted in earlier studies that such uncertainties lead to 
only a rather approximate estimation of h{x). The estimation of the scaling function associated with the bias i.e. h{x) is even 
more complicated due to the fact that the two-point quantities such as the cumulant correlators and the bias are more affected 
by finite size of the catalogs. So it is not fruitful to actually integrate the exact integral expression we have derived above 
and we will replace all line of sight integrals with its approximate values, recent study by Munshi & Jain (1999) have used 
an exactly similar approximation to simplify the one-point probability distribution function for and found good agreement 
with ray tracing simulations. We will show that our approximation reproduces the numerical results quite accurately for a 
wide range of smoothing angle. 
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Use of these approximations gives us the leading order contributions to these integrals and we can check that to this order 
we recover the factorization property of the generating function i.e. /3r,{yi,y2) = /3?)(2/i)/5»)(2/2), 



PviVuVi) = I3niyi)l3n{y2) = (2/i)/3i-h« (2/2) = T{yi)T{y2). 



(59) 



So it is clear that at this level of approximation, due to the factorization property of the cumulant correlators, the bias 
function 6,, (a;) associated with the peaks in the convergence field Ks, beyond certain threshold, obeys a similar factorization 
property too, which is exactly same as its density field counterpart. Earlier studies have established such a correspondence 
between convergence field and density field in the case of one-point probability distribution function p{5) (Munshi & Jain 

1999b), 



bn{xi)hr,{xi)bn{x2)hn{x2) = 61+4 (xi)/ii+i(a;i)6i+i(a;2)/ii+5 (2^2) • 
Where we have used the following relation between f3r, (y) and br, {x) 

bn{x)hn{x) = -— / dyT{y) exp{xy). 
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For all practical purpose we found that the differential bias as defined above is more difficult to measure from numerical 
simulations as compared to its integral counterpart where we concentrate on the bias associated with peaks above certain 
threshold. 
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It is important to notice that although the bias b{x) associated with the convergence field and the underlying density field 
is exactly same, the variance associated with the density field is very high but the projection effects in the convergence field 




Figure 3. The cumulative bias function 6(> k{9o)) is plotted against the smoothed convergence field ft((9o) = i^s- The difli'erent panels 
correspond to the different cosmological models and the smoothing angle do is fixed at 1'. The solid line at each panel correspond to 
analytical predictions from the hierarchical ansatz. The source redshift Zg is fixed at unity. 



brings down the variance in the convergence field to less than unity which indicates that we have to use the integral definition 
of bias to recover it from its generating function (see eq.(^) and eq.(^)). 

Now writing down the full two point probability distribution function for two correlated spots in terms of the convergence 
field K and its reduced version rj: 



K2)dKidK2 = PK(Kl)pK(/t2)(l + b^{Ki)^i2bK{K,2))dKldK2, (63) 

Prj{vi,V2)d-r]idr]2 = p,,(?7i)p^ (772) (1 + &^ (771)^^2^'? ('?2))rf»7irf»72 (64) 

In our earlier analysis (Munshi et al. 1999b) we found that p'^(ti) — we also noticed that ^"2 = Using these 
relations we can now write: 

b4^) = (65) 

This is one of the main result of our analysis and in the next section we will show that it is indeed a very good approximation 
to numerical simulations. 



4 COMPARISON AGAINST NUMERICAL RAY TRACING SIMULATIONS 

To compare the analytical results with numerical simulations we have smoothed the convergence field or k- map generated 
from numerical simulations using a top- hat filter of suitable smoothing angle 9o. The minimum smoothing radius we have used 
is 1' which is much larger compared to the numerical resolution length scale. In our earlier studies we found that numerical 
artifacts are negligible even for angles as small as .25'. The maximum smoothing radius we have studied is 8', which is much 
smaller compared to the box size of all models and we expect that the finite volume corrections are not significant in our 
studies. The box size is L = 166.28' for the EDS models, L =235.68' for n = 0.3 open model and L = 209.44' for D, = 0.3 
model with cosmological constant A = 0.7. The numerical outputs from more than ten realizations of these cosmological 
models were used to find the average and the scatter and were tested against the analytical predictions. For LCDM models 
we have used only one particular relaisation. 

For the determination of the bias from the convergence maps we found that the computation of cumulative bias fe(> Ks) is 
much more stable compared to the computation of the differential bias b{ns). We have therefore used the analytical predictions 
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Figure 4. Same as Figure - 3 but for smootliing angle 8 = 2'. 



1000 



100 



10 



1 

1000 



100 



10 



5CDM 



OCDM 



LCDM 




TCDM 




-0.01 

Log,„(1+K(e„)) 



Figure 5. Same as Figure - 3 but for smootliing angle 8 = 4'. 




Figure 6. Same as Figure - 3 but for smootiiing angle d = 8'. 



for &(> his) to compare against numerical simulations. The bias we have studied will represent the correlation of those spots 
in the convergence maps which are above certain threshold and is defined as: 
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The above expression relies on evaluating bias by finding out directly correlation of points in the convergence map which are 
above certain threshold and does not depend on any factorization property of bias. We have tested the functional form for 
the bias as predicted from the analytical model against numerical simulations. The analytical prediction tells us that the bias 
associated with spots where the smoothed convergence field is negative is very small and is almost zero. On the other hand 
the negative convergence spots are not biased at all, this property also holds for the bias associated with under-dense regions 
in the density field which are responsible for producing the negative spots in the convergence field. For positive peaks in the 
convergence map the bias is positive and it increases sharply as we move towards higher thresholds for Kg. Contributions to 
positive spots in the convergence maps comes from density peaks in underlying 3D mass distributions and again its remarkable 
that the bias associated with both of them are actually described by the same functional form. While we change variables from 
the reduced convergence rjs to the actual convergence field Ha we introduce a multiplicative factor of 1/kmin- As we explained 
before that this term introduces the cosmological dependence for the bias in the convergence field. This multiplicative factor is 
largest for the OCDM model and is smallest for the SCDM model and hence the bias associated with peaks in the convergence 
maps also follow the same pattern. 

The analytical results and results of numerical simulations are plotted in figures - 3-6. Comparing the analytical results 
with the results obtained from the numerical simulations we find that these predictions match very accurately. The general 
trend of increase of bias with Ks is well reproduced by the analytical models. Also we find that the analytical model is 
able to describe the relative ordering of the bias function for different cosmological models quite accurately. We have also 
computed the scatter in numerical results and found them to increase with Hs- With increasing the threshold we decrease 
the number of spots which cross that threshold thereby increasing the sample variance. The sample variance also increases 
with the smoothing length scale 9o due to decrease in the number of smoothed patches which carry completely independent 
information. To improve upon the situation we have increased the sampling by shifting the grid used to bin the convergence 
map in two perpendicular direction. Although there has been several detailed studies to quantify the effect of finite volume 
correction on the statistics of galaxy clustering, similar studies are still lacking for weak lensing statistics. However it was 
found in numerical studies clustering that the computed PDF shows a sharp cutoff at some maximum value of cell count 
TLmax which is the densest cell in the catalog. Increasing the size of the catalog will increase Umax up to which one can reliably 
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Figure 7. The analytical predictions for the bias associate with the convergence field k{9o) is plotted as a function of k((9o)- The 
smoothing angle varies from 1' to 8' as indicated. The source redshift is fixed at Zs = .5. Different curves correspond to SCDM, LCDM, 
OCDM and TCDM models. 

compute the PDF. The long exponential tail of PDF also shows large fluctuations due to the presence or absence of any rare 
dense objects in the N-body catalog. In this paper we have shown how intimately statistics of the weak lensing convergence 
field Ka and its density field counterpart are related. This will indicate that the numerical computations of the one-point PDF 
and the bias of Ks and 5 are both affected in a very similar manner due to the finite size of the catalog. It is therefore quite 
remarkable that despite the fact that we have not introduced any involved correction procedure in our measurement of bias 
or PDF the numerical results match the analytical predictions so accurately. We have not studied the factorization property 
of the bias as predicted by hierarchical ansatz, as there has not been any parallel study in this direction so far dark matter 
clustering in 3D. However we plan to extend our results to such analysis elsewhere. 

We have studied both analytically and numerically how the bias changes if we change the source redshift Zg. In our earlier 
studies of PDF we have shown that for low source redshift when one probes the highly non-linear regime directly, the PDF is 
charecterised by a long exponential tail originating from collapsed objects. Where as if we increase the source redshift we are 
adding more slabs of matter in the intervening path which are virtually uncorrelated and hence it makes the distribution more 
Gaussian. However it also increases the variance of the convergence field. We find a similar result for the bias in Ks field too, 
for low source redshift the highly dense spots which are contributing towards the nearby high KsSpots may also be physically 
be very near by and hence they themselves may be correlated more strongly than the underlying mass distribution. On the 
other hand when we increase the source redshift we increase the number of hight density spots which will produce the high 
convergence Hs spots in the map, however it is to be kept in mind that clearly not all these high density spots are in physical 
proximity in 3D space and might exist in different slabs of the underlying mass distribution. This explains why although the 
variance of Ks and correlation in Ks do increase with the source redshift but the bias associated with high Ks spots actually 
decrease with the source redshift. 

In our analytical studies we have assumed that the small angle approximation is valid and we also assumed that the 
separation angle 612 is much larger than the smoothing angle 60. This will imply that variance within the smoothed patch 
is much larger compared to the correlation between the two smoothed patches. However the separation angle should again 
be small so that we can still replace all spherical harmonics by Fourier modes. In numerical simulations we have placed the 
smoothed patches separated by an angle at least by thrice the smoothing angle. In almost all cases the separation angle is 3.5 
times larger than the smoothing angle. We have also found (as predicted by analytical results) that once the large separation 
limit is reached the bias function no- longer depends on the separation angle 812- 
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Figure 8. Same as Figure - 7 but for the redshift Zs = 2. Note that the bias associated with high k{9o) spots in the convergence field 
decreases as we increase the source redshift Zg. 

5 DISCUSSION 

The hierarchical ansatz provides a framework in which gravitational clustering is generally studied in the highly non-linear 
regime. However most of such earlier studies are directed towards understanding of galaxy clustering and were tested against 
numerical N-body simulations. However studies based on clustering properties of galaxies are more difficult to interpret as 
they always depend on a particular biasing scheme employed to relate galaxy population to the underlying mass fluctuation. 
Forthcoming weak lensing surveys will provide an unique opportunity to directly study matter distribution in the universe and 
will help us in understanding their clustering properties. The main motivation of this paper is to extend use of hierarchical 
ansatz to the study of bias as measured from weak lensing surveys. 

There have been several studies in which perturbative techniques were used to relate the convergence statistics measured 
from weak lensing surveys to the clustering statistics of the underlying mass distribution. However analytical predictions 
based on perturbative techniques normally requires a large smoothing angle. Most of early observational studies of weak 
lensing surveys will focus on smaller angular scales and hence will require a completely non- linear treatment. In this paper 
we have shown that our current level of understanding of gravitationally clustering is enough to make concrete predictions 
for two-point statistics of the convergence field. 

We have used the non-local scaling relations to evolve the matter two-point correlation function and used the hierarchical 
ansatz to relate the higher order correlation function with two-point correlation function. Combined with the generating 
function approach these two ansatze provide a powerful tool for the analysis of the weak lensing maps and its multi-point 
statistics. In accompanying papers we have computed the lower order cumulants and the cumulant correlators for various 
smoothing angle and for different source redshifts (Munshi & Jain 1999a). We have also studied the complete PDF for 
the smoothed convergence field (Munshi & Jain 1999b). In this paper we have extended these studies by testing analytical 
predictions for the bias against numerical simulations. The cumulant correlators are nothing but lower order moments of the 
bias function which we have studied in this paper. We found a very good agreement between analytical results and ray tracing 
simulations. 

The bias which we have studied here is induced by gravitational clustering. It is known from earlier studies that such 
a bicis occurs naturally both in the quasi-linear (Bernardeau 1994) and the highly non-linear regime for over-dense objects 
(Bernardeau & Schaeffer 1992, 1999, Munshi et al. 1999a,b). It is factorizable and will only depend on the scaling parameters 
X associated with collapsed objects. We have extended such analytical studies for the weak-lensing convergence field and have 



shown that similar results do indeed hold even for such cases. Our numerical study confirms such analytical claims and also 
shows that the functional form for the bias EissGciated with the convergence field is exactly same as that for over-dense spots 
in underlying mass distribution, so measurement of the bias from weak lensing observations will provide a direct probe into 
the bias associated with over-dense dark objects in the underlying mass distribution as induced by gravitational clustering. 
This will help us to separate the contribution to the bias due to gravitational clustering from the bias associated with other 
non-gravitational sources involved in galaxy formation process. 

Although our formalism (which is based on formalism developed by Bernardeau & Schaeffer (1992)) relates the bias of 
overdense cells in density field with "hot-spots" of the convergence map. Similar results can be obtained by using the extended 
Press-Schechter formalism to relate the peaks in the convergence field directly with the bias associated with collapsed halos. 
A detailed analysis will be presented elsewhere. 

The hierarchical ansatz not only provides the two-point correlation function (i.e. the bias which we have studied here) for 
the over-dense objects but it also provides a recipe for computing the higher order correlation functions and the associated 
cumulants and cumulant correlators. Our formalism which we have developed here can directly be extended to compute such 
quantities in the highly non-linear regime i.e. at small smoothing angle. However a detailed error analysis is necessary to 
study the possibility of estimating such quantities from weak lensing observations. 

Our analytical results use the top-hat filter, however it was pointed out before that the compensated filter is more suitable 
for observational studies. We plan to extend our analytical and numerical studies for other filters and results will be presented 
elsewhere. However it should be kept in mind that although nearby cells do have correlated convergence fields after smoothing 
with the top-hat window, such correlations are much less significant for the compensated filter (a property which makes it 
more suitable for the observational purposes for PDF studies). 

It is clear that the Born approximation is a necessary ingredient in the analytical computation which we have undertaken 
in this paper. The Born approximation neglects all the higher order terms in the photon propagation equation. It is then 
obvious that any statistical quantity such as the higher order cumulants and the cumulant correlators will have correction 
terms which are neglected in our approach. In the quasi-linear regime i.e. for large smoothing length scales these correction 
terms can in fact be evaluated using perturbative calculations. However such an approach is not possible to adopt in the 
highly non-linear regime, because in small angular scales where most contribution comes from length scales with variance well 
above unity the whole perturbative analysis breaks down and the only way we can check such an approximation is to compare 
analytical predictions against numerical simulations. Earlier studies of lower order cumulant correlators confirmed that such 
correction terms are indeed negligible even in small smoothing angle (Munshi & Jain 1999b). A very good match that we 
have found for the bias function in this study shows that such correction terms are indeed negligible at arbitrary order. 
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